Ratchet transport and periodic structures in parameter space 
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Ratchet models are prominent candidates to describe the transport phenomenum in nature in the 
absence of external bias. This work analyzes the parameter space of a discrete ratchet model and 
gives direct connections between chaotic domains and a family of isoperiodic stable structures with 
the ratchet current. The isoperiodic structures appear along preferred direction in the parameter 
space giving a guide to follow the current, which usually increases inside the structures but is 
independent of the corresponding period. One of such structures has the shrimp-shaped form which 
is known to be an universal structure in the parameter space of dissipative systems. Currents in 
parameter space provide a direct measure of the momentum asymmetry of the multistable and 
chaotic attractors times the size of the corresponding basin of attraction. Transport structures are 
shown to exist in the parameter space of the Langevin equation with an external oscillating force. 
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The description of the ratchet transport of particles 
in nature has become an actual and large studied prob- 
lem due to the possibility to obtain transport properties 
without external bias. To obtain ratchet transport, spa- 
tiotemporal symmetries must be broken in the system 
[1]. Ratchets have become natural candidates to describe 
transport phenomena in Brownian [2, 3j and molecular 
motors [4J , cold atoms [5j , migration of bacteria [6] , cell 
mobility in cancer metastasis [7J, granular gas [8], fluid 
transport [9j and in more general areas like classical and 
quantum physics [lOj, chemistry fill [T^] biophysics 
[13j. These are just some references in the distinct areas, 
since the actual literature related to ratchets is enormous. 

A common feature of interest in all areas of ratchets ap- 
plications is the understanding, achievement and control 
of transport. A priori, dynamical variables and param- 
eters of the system (like temperature, dissipation, noise 
intensity, external forces etc.), which control the dynam- 
ics, are deeply interconnected so that it is very hard to 
make general statements about the ratchet current (TIC) 
as a function of the parameters. The precise determi- 
nation of the nature of transport in unbounded systems 
is still not fully understood, thus it is very desirable to 
achieve and/or recognize "patterns" or "structures" in the 
parameter space which are directly connected to trans- 
port properties. Even more attractive if such transport 
structures present universal features observed in a large 
class of dynamical systems. 

This Letter analyzes the parameter space of a ratchet 
model and shows the relation between ratchet currents 
with a family of isoperiodic stable structures (ISSs) and 
chaotic domains in parameter space. In this way a re- 
markable complete connection between parameters of the 
system and the 1ZC is given, and therefore general clues 
for the origin of directed transport. To mention an ex- 
ample, one of the ISSs observed here, which has the 
shrimp-shaped form [see Fig.|3|c)], has already appeared 



in the parameter space of generic dynamical systems and 
applications. Such shrimps were found to be universal 
structures in the parameter space of dissipative systems 
like maps [T4]4T6] and continuous models [17J, among 
others. Very recently their were also observed in exper- 
iments with electronic circuits ^8j. We show here that 
such shrimp-shaped structures are also essential concern- 
ing directed transport in nature and we stress that this 
is valid for any ratchet model and applications. 

In order to show generic properties of the IZC in the pa- 
rameter space, we use a Map which presents all essential 
features regarding unbiased current [19] 

^ . f Pn+l = IPn + K[sm{Xn) + a sin(2Xn + 0)], /-^x 

where Pn is the momentum variable conjugated to x^, 
n = 1,2,...,A^ represents the discrete time and K is 
the nonlinearity parameter. The dissipation parameter 
7 reaches the overdamping limit for 7 = and the con- 
servative limit for 7 = 1. The ratchet effect appears due 
to the spatial asymmetry, which occurs with a 7^ and 
(j) ^ mn (m = 1,2,...), in addition to the time rever- 
sal asymmetry for 7 ^ 1. The TZC of the above model 
was studied [19j for fixed K = 6.5 in the dissipation in- 
terval < 7 < 1. It was shown that close to the limit 
7 = 1 the IZC arises due to the mixture of chaotic motion 
with tiny island (accelerator modes) from the conserva- 
tive case, while for smaller values of 7, chaotic and stable 
periodic motion generates the current. 

Figure [l] shows the TZC (colors) as a function of the 
dissipation parameter 7 and the nonlinearity parameter 
K. A remarkable complex structure of colors is evident, 
where each color is related to a given value of the cur- 
rent (see color bar). Black colors are related to close 
to zero currents; green, blue to purple colors are related 
to increasing positive currents while red, yellow to white 
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Figure 1. (Color online) The current (see color bar) plotted in 
the parameter space {K, 7) with a grid of 600 x 600 points, a = 
0.5, = 7r/2, 10^ initial conditions with < po >=< xq >= 
inside the unit cell (— 27r, 27r), and N = 10^ iterations. 




Figure 2. (Color online) Period-g' values in parameter space 
{K,j): green: ^ = 1, blue: q — cyan: ^ = 3, yellow: ^ = 4, 
pink: ^ = 6, red q > S and black for no period. In this 
case only one initial condition is used (xo = 0.5, po = 0.3), 
N — 10^ iterations and a grid of 600 x 600 points. 



colors related to increasing negative currents. The white 
straight line at i^T = 6.5 corresponds to the case analyzed 
recently [19j. Three main regions with distinct behaviors 
can be identified: (i) a large '^cloudy^^ background^ identi- 
fied as A in Fig.[T] mixed with black, red and green colors, 
showing a mixture of zero, small negative and positive 
currents, respectively; (ii) several structures with sharp 
borders and distinct colors, which are embedded in the 
cloudy background region and are identified in Fig. [T] as 
Bl-, Cl and Dl {Lis an integer or rational number); (iii) 
strong positive and negative currents (region E)^ with not 
well defined borders which occur close to the conserva- 
tive limit 7 = 1. At next we explain in more details these 
distinct regions by analyzing other quantities. 

Figure [2] shows the parameter space {K^ 7) for the 
period-g' from the orbits. Periodic stable motion is re- 
stricted to well defined structures while the black back- 
ground is related to the chaotic motion. This was checked 
by determining (not shown) the parameter space for the 
largest Lyapunov exponent (LE). Zero and negative LEs 
are related to the periodic motion and positive LEs to 
the black regions of Fig. [2] This already allow us to asso- 
ciate the cloudy region A from Fig. [T] with the chaotic 
motion. Thus the small portions of negative/positive 
currents (red and green clouds) are due to the chaotic 
transport, consequence of the asymmetry of the chaotic 
attractor [5j. 

Map ([1]) is periodic in x with period 27r. Thus the 
condition for a period-g' orbit is Vi — 27rm, with 

integer m = . . . , —2, —1, 0, 1, 2, . . . which represent the 
net number of left (negative) or right (positive) 27r jumps 
in X over the period. The mean momentum for a given 
period-g is then = Y^^=n~^ P'^/^ ~ ^Trm/g' which can 
be written as 27rL, so that L = m/q. Therefore L, which 
gives Pq in units of 27r, can assume fractional values (pos- 
itive or negative). All sharp borders from Bl-,Cl and Dl 
in Fig. [1] coincide with those of Fig. |2] and a direct con- 



nection between periodic stable motion and the IZC is 
given. Regarding these sharp structures we observe three 
different ISSs: Bl composed by the sequence L = 1, 2, . . . 
of dominant diagonal structures with main period q = 1 
(see green main regions in Fig. [2| which extend them- 
selves along a large range of K values. As K increases 
and 7 decreases, inside each structure Bl we observe 
1x2^ (green blue yellow ^ . . .) and 3x2^ (cyan 
pink . . .) doubling cascades bifurcations. Compar- 
ing to Fig. [1] we observe that the current inside each 
structure Bl is independent of the period q and thus, 
on the bifurcation points. At the doubling bifurcation 
cascades 1 x 2^ we observed that m 2m so that L 
(= m/q = p^/27r) remains constant inside each structure 
Bl' On the other hand, even though different isoperi- 
odic Bl structures have the same period, IZCs increase 
with L, as can be seen in Fig. [T] when 7 and L increase. 
Structures Bl are very similar to the cuspidal singulari- 
ties [20J. Since the boundary of the dominant structures 
Bl are born at period g = 1, these borders can be calcu- 
lated analytically from the eigenvalues of the Jacobian of 
the map ([T]) after one iteration. Using ^ = 7r/2, the fixed 
points from ([l]) can be calculated from p^^^ = 2iiL {L 
integer) and 27rL(7 — l)-\-K [sin {x^^"^ ) + a cos {2x^^'^ )] = 

0. The solutions are x^^^ — arctan (a^~\ ±/3^+^) (j = 

1,2) and xi^^ = arctan (a^+\ ±/3^~^) {s = 3,4) where 
a(^) = tV^ [3i^ + 8L7r(7-l)] + K and (3^^^ = 

^J^KLiril - 7) ± 2^lO[?>K + 8L7r(7 Substitut- 
ing these solutions in the Jacobian of the map ([T]) we ob- 
tain analytical expressions A(7, iC, L) for the two eigen- 
values. When A(7, iC, L) = +1 (born of period-1), we 
obtain a relation between 7, and L where period-1 
orbits are born in parameter space. The solutions for 
A(7,i^, L) — 1 = 0, for all fixed points are 7^^^ = 
1 - 3K/{S7tL) (lower border) and 7^^^ = 1 - K/{47tL) 
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(upper border). Both curves define exactly, for a given 
L, the sharp period-1 Bl boundaries of Figs. [T] and |2j 
The border of the first large dominant structure Bi is 
obtained from L = 1, the second one B2 from L = 2, and 
so on. See white lines 'y[^~^'^ , 72^""*^^ , 7i^~^^ , and 
in Figj2j Looking carefully however, all (for any V) lower 

borders do not match exactly with the simulations. 
For these regions the basin of attraction from the chaotic 
attractor is much larger compared to the basin of at- 
traction related to the fixed point, and thus the current 
related to period 1 is too small and is not observed. This 
is also exactly what is observed in Fig. 3 from [19j, where 
the left limits of the L intervals are inside the chaotic re- 
gion. As L increases, 7^^^^ 1^^^^ 1 and the 
dominant structures approach to each other more and 
more. 

The second kind of relevant ISSs, Cl, can be visual- 
ized in Fig. [1] close to i^T = 6.0 and 7 > 0.6. They also 
are ordered in a sequence of structures L = — 1, — 2, . . . 
which approach each other as 7 increases, defining a di- 
rection in the parameter space obtained by the straight 
line 7 = 0.2845i^ — 0.994925, along which negative cur- 
rents increase [see dashed line in Figjsja)]. The main pe- 
riod of each Cl is = 2 (blue) [Fig.[3p))] but, inside each 
structure, 1x2^ (blue yellow red . . .) doubling 
cascades bifurcations appears when going to the border of 
the structures, where the chaotic region is reached. For 
clarification Figs, [s]^ a) -(b) present the current and period 
as a magnification for these structure. Again we observe 
that current increases (in modulus) along the sequences 
as 7 increases. 

The last observed ISSs, Dl^ appear embedded in 
the cloudy chaotic background and present the well 
known shrimp-shaped form. For example, three con- 
nected structures (D-i) appear in red in the interval 
10 < < 12 and O.K 7 < 0.4 [see also Fig.jsfc)]. An- 
other example of shrimp-shaped structures is demarked 
in Fig. [1] by l^_i/2- These structures have a main body 
with period q = 2 and a succession of domains related to 
period-doubling route to chaos. Shrimps-shaped struc- 
tures are also distributed in sequences along preferred di- 
rection in the parameter space, as can be seen by looking 
carefully to Fig. [l] where many shrimp-shaped structures 
are hidden behind the dominant Bl structures. Such 
preferential directions in parameter space appear to be 
general properties of shrimp-shaped structures (see [14J). 
These structures are abundant in the parameter space 
as can be seen by the magnification shown in Figs, jsjc)- 
(d). As the parameter space is searched further and fur- 
ther for finer domains, a large amount of distinct ISSs 
appear, usually well organized and sometimes even con- 
nected to each other [see the connected shrimp-shaped 
ISSs in Fig.jsj^c)]. Connected shrimps have the same cur- 
rent value. Although almost all ISSs present finite T^Cs, 
some ISSs present zero currents inside. We mention one 



example which can nicely be observed in Fig.[3[a) close to 
K = 4.5 and 7 = 0.75 [compare with the cyan ISS from 
Fig. [sj^b)]. Such very interesting cases will be analyzed 
in another work. 
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Figure 3. (Color online) Magnifications of Fig. [T] (left) and 
Fig. [2] (right) showing that the ISSs are abundant inside the 
chaotic region and usually organized along specific lines and 
sometimes connected to each other. 

The region E presents larger currents close to 7 = 1 
with not well defined borders. We start by mention that 
the parameter space for the largest LE (not shown here) 
in this region is positive and thus totally chaotic. We 
also clearly see from Figs. [T] and |2] that ISSs with larger 
L start to overlap to each other when 7^1 and thus 
the number of stable periodic points in phase space in- 
creases more and more. But the interesting point is that 
apparently the periodic structures are not directly re- 
sponsible for the currents in region E. For one example 
see the preferential line along the structures C-2 • - • 
in Fig. [sj^a), where the currents increase negatively as 
7^1. If the Cl structures were responsible for the 
current in we would expect a very large negative cur- 
rent. However, this is not what happens since a current 
reversal occurs close to 7 ^ 1 (see green region in Fig. [l] 
at the end of the C_2 • . • sequence). In fact, close 
to 7 = 1 we have regions of multistability and chaotic 
motion were the dynamics, which is a mixture of a large 
number of periodic and chaotic attractors, each one with 
his own basin of attraction, is very rich and complex. 
Such multistability regions were analyzed for the kicked 
rotor in the beautiful works of ^2 The "competi- 
tion" between multistable and chaotic attractors in order 
to generate the IZC is gained by the larger basin of attrac- 
tions of chaotic attractors. Extensive numerical simula- 
tions show that close to 7 = 1 the unsharply borders do 
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not change for larger iteration times and that the basins 
of attraction of periodic orbits are very smah compared 
to the basins of chaotic attr actors. This agrees with some 
recent works [5l [H |23l [24] which suggest that the IZC in 
this region is due to the mixture of chaotic motion with 
tiny island from the conservative case. In fact, accelera- 
tor modes from 7 = 1 are responsible for the asymmetry 
of chaotic attractors, generating the currents. 

In order to show that IZCs present generic ISSs in 
the parameter space of a more general class of dynami- 
cal systems, we analyze the zero temperature Langevin 
equation: x + 7X — 5.0 [sin (x) + O.Tsin (2x — |)] — 
i^^t sin (1.0 = 0. Kt is the amplitude of the external 
time oscillating force, 7 is the viscosity and the force com- 
ing from the ratchet potential is identical from Eq. ([T]). 
Figure |4] shows the IZC in a smaller portion of the pa- 




Figure 4. (Color online) Current in the parameter space 
(i^t, X — e^"^-*) for the Langevin equation. 

rameter space {Kt^x — As in Fig. [l] larger cur- 

rents can be observed in the large blue shrimp-shaped 
ISS with main period q = 2. Again doubling bifurcations 
occur inside the ISS (not shown), as observed in Figs. |2] 
andjsj^d) for the D-i and D_i/2 ISSs. Behind the main 
shrimp-shaped ISS, a cloudy background is observed and 
also a large amount of smaller ISSs, which can be bet- 
ter resolved when the parameter space is searched for 
finer and finer domains. Small finite temperatures will 
slowly transform the sharp borders of the ISSs into un- 
sharply borders and simultaneously enlarge them since 
their attractors are more stable under noise effects than 
the chaotic attractors [25j. 

Concluding, universal ISSs in parameter space are 
shown to generate large 1ZCs and to organize themselves 
along preferential directions, which allows to make pre- 
dictions of the IZCs along the such directions. Inside a 
given structure, IZCs are independent of the period of the 
orbit. The essential property to obtain finite 1ZCs is the 
momentum asymmetry of attractors in phase space times 
the size of the corresponding basin of attraction and can 
be expressed as TIC = "^Zi^i {p)i '^i^ where Na is the num- 



ber of attractors, (p)- is the mean momentum of attrac- 
tor i and Si is the size (normalized) of attractor i. Thus 
Figs. [1] and [4] are a direct quantitative measure of 1ZCs. 
We stress that the isoperiodic stable transport structures 
Bl (cuspidal-shaped), Dl (shrimp-shaped) and Cl are 
universal patterns which should appear in the parameter 
space of any dissipative ratchet system, independent of 
its application in nature. This was ratified by showing the 
appearance of the shrimp-shaped ISSs in the parameter 
space of the Langevin equation with an external unbiased 
field. It would be very likely to observe the IZCs ISSs in 
real experiments. 
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